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Abstract: The large N A/'=4 gauge theory with quenched N=2 quark matter displays 
chiral symmetry breaking in the presence of a magnetic field. We previously studied the 
temperature and chemical potential phase structure of this theory in the grand canonical 
ensemble - here we, in addition, include the effect of an electric field which acts to counter 
chiral symmetry breaking by disassociating mesons. We compute using the gravity dual 
based on the D3/probe-D7 brane system. The theory displays two transition at one of which 
chiral symmetry is restored. At the other transition density switches on, the mesons of the 
theory become unstable and a current forms, making it a conductor-insulator transition. 
Through the temperature, electric field, chemical potential volume (at fixed magnetic field 
parallel to the electric field) these transitions can coincide or separate at critical points, 
and be first order or second order. We map out this full phase structure which provides 
varied computable examples relevant to strongly coupled gauge theories and potentially 
condensed matter systems. 
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1 Introduction and summary 

Holographic techniques [1-3] have begun to allow us to study the phase structure of strongly 
coupled gauge theories. Such studies are potentially of interest for QCD, more exotic gauge 
theories and even to condensed matter systems. Although it is somewhat hard to dial a 
particular gauge theory, we instead work in tractable models that lie close to A/" = 4 Super 
Yang-Mills theory and hope that these gauge theories show similar behaviours to realistic 
cases. 

A simple model to study is the large N J\f = 2 gauge theory made from AA = 4 
SYM plus a small number of quark hypermultiplets [4-8]. In the limit Nf <ti N quark 
loops are quenched which on the gravity side implies that one can study probe D7 (or 
D5) branes in the AdS background of N D3 branes. In the presence of a magnetic field 
(-B) the theory dynamically generates a quark condensate that breaks a chiral U{1) global 
symmetry [9]. This dynamics may be relevant to condensed matter systems [10-12] but one 
can also study it as a loose analogue of QCD-like gauge theories [13]. Of course in those 
theories the running coupling generates a scale and the chiral symmetry breaking, but the 
hard magnetic field can be thought of as the analogue of the scale Aqcd which breaks the 
conformal symmetry and "allows" the strong dynamics to generate a condensate (c). 

In [13] we computed the phase diagram for the theory with massless quarks, in a 
magnetic field, in the temperature chemical potential plane (T-fx plane) which we show in 
Fig 1(a). It has considerable structure. Whilst the magnetic field favours chiral symmetry 
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breaking the temperature [14] and chemical potential [15] favour the chirally symmetric 
phase. This tension leads to interesting phase structure. There are three phases - the 
high T, fj. phase with the symmetry restored (xS), finite density (d ^ 0) and no stable 
mesonic states; an intermediate phase in which chiral symmetry is broken (xSB), finite 
density {d ^ 0), but mesons are still unstable; and finally a low T, /u phase with broken 
chiral symmetry (^SB), no density {d = 0) and stable mesons. On the gravity side these 
phases are described by a flat D7 embedding, an embedding that spikes on to the black 
hole horizon, and a solution that lies outside the horizon (See the insets in Fig 1(a) for 
schematic plots.). These phases are linked by transition lines that are in places first order 
(blue lines) and in places second order (red lines) - critical points link these behaviours. 
Time dependent aspects of these transitions have also been studied in [16-19]. 

In this paper we wish to study how robust the phase structure is to changes of pa- 
rameters. In particular we will introduce the additional parameter of an electric field (E) 
parallel to the magnetic field. Note this case is simpler than when the E and B fields are 
perpendicular because no Hall current forms. We do allow for the induced current in the 
direction of the electric field. The effect on the theory of an electric field has been stud- 
ied previously for probe branes in [20-26]. Although the electric field continuously acts 
on the quarks an equilibrium configuration is nevertheless reached where energy is being 
dissipated to the bulk A/" = 4 plasma. The quarks and anti-quarks have opposite charge 
when interacting with the electric field so the field tends to loosen the binding in mesonic 
bound states and, if it is strong enough, to disassociate the mesons. With an electric field 
present in the theory a singular shell develops in the gravity description. If the D7 brane 
passes through this shell its action becomes imaginary and in order to keep the action real 
one must turn on the appropriate electric current (J). The singular shell plays the role of 
an effective horizon for world volume meson fiuctuations [27]. If the probe brane touches 
the singular shell one expects that fiuctuations of the brane there must be in-falling and 
the spectrum will resemble a quasi-normal mode spectrum describing mesons that have a 
complex mass. In other words the electric field has acted to disassociate the meson in a 
similar fashion to how temperature melts the mesons [28]. 

The field theory phase structure is determined holographically by comparison of the 
classical bulk field configurations. In this paper, we introduce three bulk fields L{p) (em- 
bedding scalar), At{p) (gauge field), and Ax{p) (gauge field) in the D7 brane world volume 
with given background parameters {E,B,T). By fixing the asymptotic values of fields in 
the UV (large p) as (L — t- m (quark mass). At ^ p (chemical potential), A^ — ?• 0), we 
look for the sub- leading behaviors of fields at large p: c (condensate), d (density), and 
J (current), which are determined by the bulk DBI dynamics. Therefore, our problem is 
classifying the phases by three quantities (c, d, J) in the 5D space (T, B, E, m, p). It turns 
out that, because of a scaling symmetry, we can scale all variables by B, which reduces 
our phase space to 4D. Since we are interested in spontaneous chiral symmetry breaking, 
we will choose m = 0. Our phase space becomes 3D (T, /i, E) and we will classify this 
space by 8 possible states consisting of the three order parameters (c, d, J) being "on or 
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(a) T-/X phase diagram. The inset diagrams are (b) T-E phase diagram. The inset diagrams are 

representative probe brane embeddings (dotted representative probe brane embeddings (dotted 

lines), where a black disk represents a black hole. lines), where a red arc represents a singular shell. 




(c) fi-E phase diagram. The ^j.{E) axis agrees with (d) E-^-T phase diagram. The curves show the 

the /i(i?) axis of Fig 1(a) (1(b)) and are shown by phase boundaries at fixed T {T — 0, 0.15, 0.21, 

the same color. Thus embeddings are the same as 0.232, 0.24 from the bottom). See Fig 7 for more 

1(a) and 1(b) without the black hole. detaU. 



Figure 1. The phase diagrams of the massless JV ~ 2 gauge theory with a magnetic field. First 
order transitions are shown in blue, second order transitions in red. The temperature is controlled 
by the parameter T, chemical potential by fi and electric field by E. 

off" . Among them, only in the c^O,d=J = phase are stable mesons allowed^ . T, fi, E 
tend to turn off c (c = 0) and turn on d,J {d,J ^ 0), and so oppose B. Due to these 
competitions between E, T, /i and B a rich phase structure is constructed. 

One of our main results here is that three phases (corresponding to three embedding 
types) are also present in the T, E plane as we summarize in Fig 1(b). With a large (small) 
electric field, the theory becomes a conductor (insulator) and is in xS (xSB) phase. At an 
intermediate electric field, the system is a chiral symmetry broken conductor. Note that 
the orders of the phase transitions and the positions and presence of the critical points 

^We should caution that when studying phase diagrams one is always limited by the states allowed in 
the analysis. We do not study the effects of the parameters on the squark potential for example which is 
likely unstable in the presence of a chemical potential [29] . The electric field could also potentially generate 
phenomena beyond chiral symmetry restoration, density creation, current induction, and meson melting 
but our results should stand as a starting point for exploring such extra phases should they exist. 
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vary relative to the T, /i case (Fig 1(a)). 

Interestingly, the ^-E phase diagram (Fig 1(c)) shows a very different structure from 
the T-fi or the T-E phase diagram. At zero temperature and finite E and B field, the 
contribution to the action from the density is canceled by the contribution from the induced 
current, which is a function of density. Consequently the free energy of the system is 
independent of density, which was pointed out in [22] in the zero B field case. Thus the 
system is essentially a zero density system and the phase diagram is independent of the 
chemical potential. At first this may seem at odds with Fig 1(a) since the n axis has 
structure present. In fact there is a first order transition at T = between the fj, axis and 
the rest of the E,fi plane. There are two limits approaching the ^u-axis: (1) E = and 
then T — )• (2) T = and then E ^ 0. These two limits are different and only the former 
is a continuous limit to the /i-axis. 

Indeed by computing the phase structure through the entire T, fi, E volume (e.g. com- 
puting the fi-E diagram at various fixed T) we show the smooth evolution through that 
volume connecting the three surface planes. There are interesting movements of the critical 
points and phase boundaries. These results are summarized in Fig 1(d) which we discuss 
in much more detail in section 6. We also see the development of the first order transition 
between the E-fi plane and the //-axis as T — t- 0. Two of the missing states from Fig 
l(a)-l(c), (c 7^ 0, d / 0, J / 0) and (c = 0, (i 7^ 0, J / 0) from among the 8 possible states, 
are found in the 3D bulk of Fig 1(d). 

These models with their varied behaviours in the B, T, /x, E volume can hopefully serve 
as exemplars, or templates from which to find exemplars, for different phase structures in 
physical theories. 



2 The holographic description 

The AA=4 gauge theory at finite temperature has a holographic description in terms of an 
AdSs black hole geometry (with N D3 branes at its core) [1-3]. The geometry is 

ds^ = ^{-fdf + dx2) + 4^dr2 ^ ^2^5^2 ^ (2.1) 
where = iirgsNa"^ and 



r 



4 

/ := 1 - ^ , rH:= ttR^Tft • (2.2) 



r4 



Here rn is the position of the black hole horizon which is linearly related to the dual field 
theory temperature Tpx- 

We will find it useful to make the coordinate transformation 



y/2w:= \ r"^ + J r'^ - r% , (2.3) 



j.2j ^2 

with = th- This change makes the presence of a flat 6-plane perpendicular to the 
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horizon manifest. We will then write the coordinates in that plane as p and L according to 



w = yfP^+I? , p := u;sin6' . L := wcosO , (2.4) 

The metric is then 

2 7->2 

ds^ = -^{-gtdt^ + gxdx^) + -^{dp^ + p^d^l + dL^ + L^dnl) , (2.5) 

where 

9t ■■= 4. 4^ 4 N ' 9x ■■= 4 • 2.6 

Quenched {Nf <^ N) J\f=2 quark superfields can be included in the J\f=4: gauge theory 
through probe D7 branes in the geometry [4-7]. The D3-D7 strings are the quarks. D7-D7 
strings holographically describe mesonic operators and their sources. The D7 probe can 
be described by its DBI action 

Sdbi = -Td7 j d'^v'-det(P[G],6 + 2-Ka'Fab) , (2.7) 

where P[G\ab is the pullback of the metric and Fab is the gauge field living on the D7 
world volume. We will use Fab to introduce a constant magnetic field (eg F12 = —F21 = 
B/{2Tra')) [9], a chemical potential associated with baryon number At{p) 7^ [15, 30] 
and our crucial extra ingredient here an electric field parallel to the magnetic field (Fos = 
~-^30 = -E/(27ra')) [20-23]. We will also allow for the possibility that the electric field 
induces a current in the z-direction by including A^- 

We embed the D7 brane in the t, x, p and directions of the metric but to allow all 
possible embeddings must include a profile L{p) at constant Oi. The full DBI action we 
will consider is then 

S = J deC{p) = €3 J dtdx^ J dp C{p) , (2.8) 
where £3 is a volume form on the 3-sphere. Here 

£:=-A'p=fl-=iV[7(V^|VT^l (2.9) 





and A/" := NfT^-j. At large p, for fixed E and B, the fields behave as 

c d J 

L~m + ^ + ..., + ^ + (2.10) 

pZ pZ pZ 
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where m is the quark mass, c the quark condensate, ^ the chemical potential, d the quark 
density and J the current. Since the action is independent of At and A^, there are conserved 
quantities d and J ^:= ;5^^ • These relations can be inverted to express A'f. and 

A' in terms of d and J as 



27ra'A', 



d w"^ — 



T 4 I 



27ra'A/' 



w 



H 



Q:-- 



1 - ^2_R4 



2 «;4(u,4+^4^) / 



V" 



This is used to express the Legendre transformed action in terms of d and J: 



Slt = S 



dfFpt 



6S 



SF, 



pt 



6S 



j j dtdx^ j dp Clt{p) , 



where 



K 



1 - 



•2 c>4,„4 



VK{1 + L'2) 



+- 



4 , 4 \ 2 



d^ 



4 n2 



/O + 



Note that the first factor of K changes sign at Ws, 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



ER^ 



+ 



(2.15) 



which defines a singular shell with a radius li;^. At zero temperature {wh = 0) the singular 
shell forms at ^/ER, and at zero E, it merges into the horizon Wg = wh- Note that the 
singular shell does not depend on density. In order to make the action regular, the second 
term of K should change sign at the singular shell. This condition determines the current 
J and conductivity a by Ohm's law: 




(2.16) 



d^ 



+ 



R* 



where ps is the p coordinate where an embedding touches the singular shell. This is still 
a function of quark mass m after all other parameters are fixed. Inspite of the way we 



-6- 



write (2.16) the current is non linear in E since Ws and ps are functions of E. a has two 
contributions. The first term is from net charge carrier density, d, and the second term 
is from pair produced virtual charges. Interestingly, the conductivity by pair produced 
charges are enhanced by B. The more general conductivity for arbitrarily angled constant 
E and B was obtained in [21, 26] in a different coordinate system^. By plugging (2.16) 
into (2.13) and rescaling we have a dimensionless Lagrangian Clt'- 



K := il 



WbW 

E^w^ 

(y;4 _ 2-4)2 



K(l+L'2) 



(2.17) 



w 



{w^-T*y{w^ + T^)wl {w^-T^y{T^ + wt)\ 

where we rescaled 

{u},L,p) ^ rVb {uj,L,p) , {d,J)^{RVBfM2TTa{d,J), E ^ BE , 



,(2.18) 



(2.19) 



and define T = wh for notational clarity. The Lagrangian Clt will be our starting point 
for the numerical analysis in the following sections. 

The chemical potential is obtained by integrating A[ (2.11) from the horizon to the 
boundary 



d 



w 



rp4 



PH 



u;4 + r4 



l-E^^^;j^] (l+L'2) 



\ " (U.4+T'4) 



j2 w^jw^+T^) 



+ 1 + 



^4 



(2.20) 



p- 



where At{pH) = 0. 

In the following sections we will present our results on various aspects of the phase 
structure of this theory. Until the final section we will concentrate on the case of massless 
quarks where the U (1) symmetry in the dVti direction is a good UV symmetry of the 
theory. Also here and below we will express all our dimensionful parameters in units of 
the magnetic field B to the appropriate power (see (2.19)) - in other words we will use the 
magnetic field as the intrinsic scale in the theory. 



3 B, T, jji Phase Diagram 

In this section we review previous results [13] on the theory without an electric field present. 
The influence of B, p and T on the D7 embedding can be understood qualitatively as 
follows. With none of these terms present the embedding of the N = 2 theory for massless 
quarks is just flat, L(/9) = - c is zero and there is no quark condensate generated to break 
the U{1) symmetry. 

^The conductivity of other models have been obtained by the same method. See for example [31] and 
references therein. 
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(a) The D7 brane embeddings (b) m- c diagrams (c) the Free energies 

Figure 2. Exemplary plots in the presence of a magnetic field at finite chemical potential and 
with temperature T = 0.2156. (Parameters are scaled by B.) 



Temperature is represented by an attractive black hole at small p and also favours a 
flat embedding. The presence of density d leads to solutions that spike to the origin (or 
on to the black hole horizon with temperature present) - the spike can be thought of as 
strings linking the D3 and D7 branes and represent the density of quarks explicitly. Note 
that at finite chemical potential there can also be solutions with At = p. and d = - the 
action for such solutions is independent of /.f and hence so are the embeddings. 

Non-trivial embeddings are generated by the B field. We can see from (2.9) that the 
last square root term grows if the embedding approaches w = Q when there is a S field. 
The embeddings tend to curve off axis to avoid the origin leading to a non-zero quark 
condensate c and spontaneous breaking of the U{1) symmetry. 

The interesting phase structure of the theory is generated by the competition between 
the B field wanting to break the U{\) symmetry and d and T which seek to restore it. As an 
example we show in Fig 2(a) solutions of the Euler-lagrange equations for the embedding 
for a case with ^/B > T > fi - the red embedding is flat and preserves the symmetry, the 
yellow embedding spikes on to the horizon and the dark blue embedding breaks the U{1) 
symmetry and has d = 0. The full phase diagram is generated by competition between 
these three types of embeddings. 

To decide which of the embeddings represents the true vacuum one may compute the 
free energy = —S evaluated on the solutions. This plot is shown on the right hand side 
of Fig 2(c). Here we see that we have chosen a transition point where the red and dark blue 
curves are degenerate - for larger T the flat embedding is preferred. Here this transition is 
associated with chiral symmetry restoration (c changes from non-zero to zero) and melting 
of the mesonic degrees of freedom (the brane moves from off the horizon to on it). It is a 
first order transition. 

The transition's presence and position in parameter space can also be seen from the 
middle plot in Fig 2(b). Here we plot the condensate c against the quark mass m. To 
obtain this plot we must find solutions for the embeddings that lie off axis as p — >• oo such 
as the green and orange curves shown. We can see the presence of the three solutions 
discussed at m = in the middle plot. To see there is a phase transition at this point we 
can use a Maxwell construction. The quark mass and condensate are conjugate variables 
and therefore the area between the c axis and the curve represents the free energy difference 



-8- 



between two points on the curve. At this point the two segments of the curve generate 
equal area and the first order transition between the red and blue points is predicted. 

In [13] we fleshed out in this way the full phase diagram of the theory which we show 
in Fig 1(a). We quickly present the results here since the analysis is a subset of the work 
in section 6 below. In the bottom left segment the symmetry breaking, black hole avoiding 
embedding is the vacuum; in the central section a spike embedding is preferred; and to 
the right the symmetry preserving flat embedding is the vacuum. The low T, // phase has 
chiral symmetry breaking, no density and stable meson states; the central regime chiral 
symmetry breaking, finite density and melted mesons; in the high T, fi phase the symmetry 
is restored and the mesons melted at finite density. The order of the transitions is marked 
as well and each transition has periods where it is first and second order, linked by critical 
points. 



4 B, T, E Phase Diagram 

The rich structure of the B, fi,T phase diagram leads one to ask how generic it is? The main 
goal of this paper is to introduce an additional parameter that favours chiral symmetry 
preservation to see how sensitive the phase diagram is to a change of parameter. We will 
use electric field (E) as that new parameter. 

4.1 B, E at zero temperature 

As a first example lets consider the system with E and B but no T ox ji. The Legendre 
transformed Lagrangian is 




(4.1) 



As has been discussed in (2.15) and (2.16), there is a singular shell at Wg = V E and the 
current is given by 

J ^ ^V^+wM ^ ^(i + ^2)^cos6 Os . (4.2) 



One can therefore find the embeddings that end on the singular shell by shooting out 
(in) from the singular shell with J determined by the point on the shell one is shooting 
from. We then seek amongst such solutions for those that satisfy L — )■ (m — )• 0) as 
p — )• oo to find the massless (mass m) embeddings. Generically we again find three types 
of solutions - (1) embeddings that curve off axis and reach p = missing the singular shell 
(when E <^ B); (2) embeddings that curve off axis and pass through the singular shell 
(when E ^ B); (3) the flat embedding L = 0. The schematic plots of these three cases are 
shown in the inset of Fig 1(b), where the black disk should be ignored at zero T. In Fig 
3(a) we show some sample numerical embeddings ending on (2) and off (1) the singular 
shell. 

All non-flat embeddings that pass through the singular shell have a conical singularity 
at /3 = 0, whose precise interpretation is unclear and discussed in [22, 23]. The conical 
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(a) Embeddings ending on and off the singular shell (£ = 0.1, m = 0). 
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(b) The condensate, c vs mass, m for given E values. Curves shift to left as E decreases. 
The red (blue) part corresponds to the embedding touching (missing) the singular shell. The 
inset on the right is a zoom-in around the first order phase transition. The colored points 
correspond to the colored embeddings in Fig 3(a). 



Figure 3. Embeddings and c-m diagram at finite E parallel to finite B, but no density and 
temperature. 



singularity is most likely a reflection of the energy being injected by the electric field being 
sunk into the gauge background through stringy physics representing the quark interactions 
with the = 4 Yang-Mills fields. 

The Fig 3(b) shows the condensate, c, vs mass, m, plot for given Es. One point (c,m) 
in the plot corresponds to one embedding since it gives a complete initial condition for the 
embedding equation (a second order differential equation). For example, the colored points 
in the inset of Fig 3(b) correspond to the embedding in Fig 3(a) of the same color. 

From Fig 3(b) we can determine the phase structure as follows. 

For a larger E (relative to B), the c-m curves tend to be pushed to the right. So if we 
focus on m = case, the only available point at large E is (m, c) = (0, 0), which corresponds 
to the fiat embedding. The fiat embedding preserves the U{1) chiral symmetry, so the 
system is chiral symmetric (xS). Since the flat embedding necessarily crosses the singular 
shell, there is no stable meson but there is a current (4.2) with (9^ = 0: J = y'{l + E'^)E, 
which is the maximal current for a given E. i.e. The system is a conductor. 

As E is lowered the c-m curve shifts to the left. At E' ~ 0.2054 a new solution for m = 
appears, whose c is finite and the ground state corresponds to this new solution. This is a 
second order transition because the second solution for c at m = moves smoothly away 
from the c = embedding. Because of the finite c, the embedding is curved and breaks 
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the U (1) symmetry, so the system is in a chiral symmetry broken (xSB) phase. However 
the embedding still passes through the singular shell, so it is a conductor with no stable 
meson. 

Finally, around E 0.1 the small 'S' shape structure connecting the red part (singular 
shell touching embedding) and blue part (a singular shell missing Minkowski embedding), 
meets the m = line, which is zoomed in, in the inset of Fig 3(b). It shows a typical 
first order phase transition structure and by a Maxwell construction we can pin down the 
transition point. So as E decreases, the green point should move to the gray point with 
a discontinuous condensate jump (see the zoomed in inset). The blue part (or the gray 
points in the inset) corresponds to the Minkowski embedding and the system is a xSB 
insulator with stable mesons. 

One would like to match this picture to a computation of the free energy. Naively it 
seems one should just compute the original action (before Legendre transforming) evaluated 
on these solutions. However, there is a subtle point related to the conical singularity^ of 
the embedding at p = and also a log divergence of [20, 24, 26] at the horizon. Thus we 
need to add some boundary term to take care of the singularity at /o = or at the horizon. 
This boundary term will also contribute to the free energy so must be taken into account. 

However, this surface term does not change the equation of motion. Furthermore, 
as far as the embedding dynamics is concerned, the singular shell position has the same 
singular structure as a black hole horizon and the embedding outside a singular shell is 
independent of the ones inside the shell. So the c-m plots, based on the classical embedding 
outside the shell, are valid regardless of the additional boundary terms at the IR boundary. 

These solutions correctly show us the maxima and minima of the free energy as a 
function of E. The discussion above is the only consistent picture with the c-m plots so 
we can be confident of its validity. For this reason in what follows we will focus on the 
c-m plots (and also when chemical potential is present we will track the quark density) to 
determine the phase structure. A similar philosophy was used in [22, 23]. 

However, it would be interesting to identify the correct boundary term and compute 
the consistent free energy graph for our c-m plot (we have not been able to so far). To 
identify it, in principle, one should start with the time-dependent back-reacted system, 
since the singularities are related to time-dependent energy loss of the system and its 
effect on the background adjoint matters. However one may also be able to introduce 
an "effective" boundary term. Our c-m diagram would be a good guide to figure out 
the correct boundary term and free energy or could even be used as a rule to determine 
it, since sometimes thermodynamic consistency plays a complementary role in AdS/CFT 
applications. 

^If we consider the finite density system this conical singularity disappears. However another singularity 
seems to appear for At as discussed in section 5. 
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(a) T = 0.24. E = 0.0001, 0.049, 0.1 from left. (b) T = 0.30. E = 0.0001, 0.100, 0.180 from left. 
Figure 4. c-m diagram at finite E parallel to finite B at finite (high) temperature. 



4.2 B, E at finite temperature 

We can now extend our analysis to include temperature straight forwardly. From (2.17) 
with d = 



K 



K(l + L'2), 



(y;4 _ 2-4)2 



4 \ 2 



(^4 _ ^4)2 



(4.3) 



and 




J = E^ 



/^4 + (^4+^4)2 

Ps 



(4.4) 



It is apparent that even with non-zero T there remains a singular shell - it always lies 
outside T for any T (2.15). Requiring regularity of the onshell action allows us to fix the 
current J (2.16). If the embedding does not touch the singular shell ps = 0, then there is 
no current. 

We shoot out to obtain the embeddings as a function of E and T at fixed B. The 
process is laborious - we plot the evolution of the c-m plot on fixed T trajectories as we did 
at T=0 in the previous section. There are three types of c-m plot. At low temperature, it 
is similar to Fig 3(b). At high temperature two qualitatively different structures appear as 
shown in Fig 4. As the temperature increases, the T effect dominates the E effect, which 
is visualized in the c-m diagram as follows. The curve near (m, c) = (0, 0) is curved to the 
left (Fig 3(b)) at low T but curved to the right (Fig 4(a)) at high T (similarly to Fig 2), 
showing competition between E and T. Both are an attractive effect from the embedding 
dynamic's point of view, but the T driven 1st order attraction is so strong that the E driven 
2nd order smooth attractive effect cannot be realized. At very high T the repulsive effect 
of B is completely suppressed and the only allowed embedding is a flat one (Fig 4(b)). 

The resulting phase diagram is Fig 1(b). For T < 0.233, three regions, as at zero T, 
exist: as E increases, the phases changes from a xSB insulator (stable meson) phase to a 
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xSB conductor at a first order transition. There is then a second order transition to a xS 
conductor. Above T ~ 0.233 the intermediate region xSB and conductor phase disappears 
and the chiral symmetry restoration and insulator-conductor transition happen at the same 
time. It is a first order transition and exists in the temperature range 0.233 ^ T < 0.25. 
Finally for higher temperatures T > 0.25, the system becomes a chiral symmetric conductor 
for any finite E. 

There are distinct features of the phase diagram though from the T, ji phase diagram 
in Fig 1(a). The insulator-conductor transition is first order along its whole length. At 
finite mass and zero magnetic field, this insulator-conductor transition was also shown to 
be first order in [22, 23]. The chiral symmetry restoration phase transition is second order 
along all its length from the critical point where it joins the insulator-conductor transition. 
The ability to reproduce different phase structures is interesting and potentially useful if 
one wanted to use these models as effective descriptions of more complex gauge theories or 
condensed matter systems. 



5 B, ji, E at zero temperature 

We now turn to the E, plane at fixed B and zero T where the behaviour appears somewhat 
different from those of the planes so far discussed. 
From (2.17) with T = 



ZiT = ^/(^TX^./(l-^ 



(5.1) 



and 



(5.2) 



where Wg = v E. Note that this formula for J is only valid for finite E, because the current 
J is introduced to make a sign change when ^1 — changes sign. If ii^ = we wouldn't 
have any reason to introduce J. 

If we now substitute J back into Clt then the d dependence explicitly vanishes, as 
also observed in [22] for the zero B case. 



p6 + -^-(l + ^2)^COs6 



It is the same action as (4.1), but the physics could still be different because at finite density 
the boundary condition for the embedding is different from zero density. Furthermore the 
current (5.2) looks different from the d = case (4.2) through its d dependence (explicitly 
and implicitly through Og). However, at zero temperature and finite E^ it turns out that 
the density is always zero, which can be shown as follows. From (2.11) or (2.20), at zero 
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Figure 5. d vs T for four points {E, fi) = (3, 10), (0.3, 10), (3, 1), (0.3, 1). d ^ as T ^ 0. 



A', = d 



\ 



-d? 



1 + 



(5.3) 



where w = \J L"^ ^ and the square bracket is the current which must be zero at = 0. 
So, for non-zero the density cancels in the denominator. Let us first consider a fixed 
non-zero E and the flat embedding and take the hmit p — )• 0. 



A' ^ ^ 



(5.4) 



Thus fi d X a d-independent integral that diverges near p = 0. There is, therefore, no 
way to get a finite /i from a finite d. The only available density is exactly zero and if density 
is zero we shouldn't use the relation (5.3). Any constant At is an available solution so any 
constant /i is allowed. 

To confirm this analysis, we numerically evaluate density at small temperatures for 
four sample points, {E,p) = (3, 10), (0.3, 10), (3, 1), (0.3, 1), in Fig 5. The density indeed 
vanishes as T goes to zero. Turning on larger chemical potential does not change the 
tendency: a 10 times larger chemical potential vanishes 10 times faster (compare ;U = 10 
and p = 1). Electric field does not affect this much: E = 3 and E = 0.3 at p = 1 are 
indistinguishable in Fig 5. Thus we find that vanishing density at T = is consistent with 
the limit T 0. 

Given this argument it is worth checking how nontrivial results come from the same 
expression on the p axis in Fig 1(a), where T = E = 0. Let us consider E = not E ^ 0. 



A[ = d. 



d^w"^ + (1 + w^)p^ ■ 



(5.5) 
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For a flat embedding and near p = 



4-1-^, (5.6) 

We learn that ;U ~ d x a d-dependent regular integral. So here there is a non-trivial 
/i, d relation. We should be careful with the d — )• limit. Then the d-dependent integral 
diverges as d gets smaller. However it turns out that it is less divergent than 1/d so we 
get // = as expected. If we consider the spiky embedding then the singular integral will 
fall as ~ 1/d due to the divergence in L'. Here we can have a nonzero chemical potential 
when d — 7- 0. This occurs at the first transition point (from the solution lying outide the 
singular shell to the spiky embedding) shown on the fi axis (zero T) in fig 1(a). 

Notice that the T = E = theory and the T = 0, S limit are distinct. The 
T = 0, E,iJ, plane for any finite E has zero density. However, the strict E = axis does 
have density present for the spike and flat embeddings - we show this with the yellow and 
green lines respectively on the E = axis in Fig 1(c). Since there is a jump in the density 
off the E = axis at T = there is formally a first order transition with increasing E, 
which is expressed by the blue line near ii^ = in Fig 1(c). (note the transitions on the 
E = axis are second order though - two red dots in Fig 1(c)). In the next section we will 
approach the T = plane from positive T to confirm this picture. 

The surprising aspect of this result is that, at zero T, and with even an infinitessimal E 
present, density is not generated no matter how large the chemical potential /i is. Therefore 
the flat, chirally symmetric configuration is not favoured for very large /i when a small E 
is present. 

These conclusions are certainly correct within the DBI analysis presented here. If the 
reader wishes to seek additional physics that might generate density at T = 0, E ^ and 
more simply connect the phases at zero and infinitessimal E, then one might be able to do 
that through additional boundary terms at the origin (this is where singular behaviour also 
enters in the //, d relation). Presumably at this point the physics associated with the sink 
of the energy being injected by the E field should be better understood. Equally the DBI 
action may not be valid for this system due to a divergent gauge field near the origin [27]. 
Resolving this issue is beyond our DBI analysis here. 



6 The Full B, T, E, /i Phase Structure 

Our final task is to complete the phase structure analysis by extending it to the full E, /i, T 
volume at fixed B. The number of embeddings that must be analyzed on any fixed plane 
through this space is already large so we restrict ourselves to looking at some representative 
slices that will be sufficient to reveal the structures present. 

In particular we will study fixed T slices and draw the phase structure in the E, fi 
plane. The embedding equations are now given by the full forms of (2.17). To determine 
the presence and nature of a transition it is sufficient to track any operator of the theory. 
We have found it easiest on these planes to plot the density d against the chemical potential 
i.e. we use d as our order parameter. 
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(a) 2nd — 1st order 
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(b) 1st — >• 2nd order 



T = 0.2.12, E = 0.04 



(c) 2nd — i- 2nd order 



T = 0.200. E = 0.085 
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(d) 2nd 1st — >• 2nd order 



(e) 1st order 



(f) 2nd order 



Figure 6. Representative n-d plots. The black line on the fj, axis corresponds to a embedding that 
misses the singular shell; the orange curve to a spike embedding that ends on the singular shell; 
and the red curve to a flat embedding. The transition points are shown by the vertical dotted lines. 
Subcaptions are the order of transition as ^ increases. 



In Fig 6 we show sample plots from which each of a first and second order transition 
can be identified. There are in total six types of d-fi plots. In each of the figures the black 
line on the /i axis corresponds to a chiral symmetry breaking embedding that misses the 
singular shell; the orange curve to a spike embedding that ends on the singular shell; and 
the red curve to a flat embedding. It can be seen from the plots whether each transition 
between embeddings is smooth and hence second order, or whether there is an 'S' shaped 
structure so that one expects a first order transition. The transition points are shown 
by the vertical dotted lines. We have used these techniques on constant E lines on each 
constant T plane to determine the transition places and orders. Fig 7 is constructed from 
sequences of d-fi plots (Fig 6). For example, Fig 7 shows the following evolution of the d-fi 
plot as E increases. 



Fig 7 
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In Fig 7 we show six slices through the volume at varying T. Starting at high tem- 
perature the theory lives in the chirally symmetric phase with unstable mesons and the 
material is a conductor. As the temperature falls to T ~ 0.25 the first order transition to 
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(a) T = 0.232 (b) T = 0.224 (c) T = 0.222 




(d) r = 0.210 (e) T = 0.200 (f) T = 0.003 (Notice that the 

scale of ^ is 100 times larger) 



Figure 7. Phase diagrams in the E plane at various vahies of temprature showing the phase 
structure. The solid blue lines are first order transitions, the red lines are second order. 

the chiral symmetry breaking, stable meson, insulator regime begins to appear in the ^-E 
plane around = E = 0. That transition then expands away from the origin and remains 
briefly first order. 

Our plot for T = 0.232 shows the first interesting structure. Two areas in the plane 
grow out from the first order line bordered by additional second order transitions. In these 
areas the theory is in a chiral symmetry breaking but conducting phase. The critical points 
where the first and second order transitions meet migrate inwards along the first order 
boundary from each axis as the temperature falls and eventually pass each other as shown 
for T = 0.224. At temperatures below that point there are constant E trajectories across 
which there are three transitions. A second order transition from conductor to insulator, 
a first order transition between two spike embeddings and finally a second order chiral 
symmetry restoration transition. An example of the relevant density chemical potential 
plot for this case is in Fig 6(d). 

Our plot for T = 0.222 shows the next key transition. The first order line between the 
two spike embeddings breaks making the spike embedding phase continuously connected 
although there is a remnant of the first order transition ending free at a critical point. 

The first order transition near the E = phase boundary then diminish as temperature 
is further reduced retreating towards the axes - see the plot for T = 0.21. It has totally 
disappeared by T = 0.2. However the first order transition element on the conductor 
insulator transition grows from the E axis as temperature decrease. 

The final interesting feature begins to appear in the plot at T = 0.21 where the phase 
boundaries have begun to deform. They expands out to large fi very rapidly. Note that 
the scale of the fi axis at T = 0.03 is 100 times larger than the others. If we drew it on 
the same scale as the others it would look like Fig 1(c): the /i-independence of the T = 
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limit is starting to be seen. At zero T (Fig 1(c)) there is a first order transition between 
the fji axis and the rest of the /i, E plane - here we see that that forms as the second order 
boundaries are pressed onto the /i axis. 

These results have been incorporated into the 3D plot in Fig 1(d) which summarizes 
the full and rich phase structure. 

7 Quark Mass 

The analysis above has been purely for zero quark mass. We have not performed an 
analysis of the introduction of a quark mass, however, we note here that an immediate 
consequence of introducing a small quark mass is that the second order chiral symmetry 
restoration transitions (the outer red lines in our figures) become crossovers. We would 
expect some remnant of the first order segment of this transition to remain associated with 
a transition at which there is a discontiunuity in the quark condensate even in the infinite 
mass limit where the theory becomes the massive N = 2 theory (we observed this in the 
E = Q limit in [13]). The insulator conductor transition remains with quark mass and will 
again become that of the N = 2 theory at large mass. At finite mass and temperature 
but without magnetic field, the insulator-conductor transition was shown to be first order 
in [22, 23]. 

8 Summary 

We have explored the phase structure of the N = 2 gauge theory whose dual is the D3/D7 
system. A magnetic field, tries to induce chiral symmetry breaking. An electric field, 
-E, tries to disassociate the mesons of the theory and makes it a conductor. Finite density, 
d, (or chemical potential /i) and temperature, T, each favour melting of the mesons. The 
competition between these effects lead to a rich phase structure. 

In Fig 1(a) we display the T, phase plane for the massless theory at fixed magnetic 
field. There are three phases - at low T, fj, a chiral symmetry breaking phase with stable 
mesons; at intermediate values a chiral symmetry breaking phase with unstable mesons; 
and at large T, fj, a chirally symmetric phase with unstable mesons. The transitions between 
these are a mix of first and second order transitions linked at two critical points. 

In Fig 1(b) we show the T, E phase plane for the massless theory at fixed magnetic field. 
There are again three phases - at low T,E a chiral symmetry breaking phase with stable 
mesons which acts as an insulator in the presence of a small electric field; at intermediate 
values a chiral symmetry breaking phase with unstable mesons which is a conductor and 
sustains a current in the presence of an electric field; and at high T,E a conducting chirally 
symmetric phase with unstable mesons. The transitions between these are again a mix of 
first and second order transitions linked at one critical point. 

The E,iJ, phase plane has rather different structure (Fig 1(c)). In particular for any 
finite E the plane is independent and density is zero. At low E we have a chiral symmetry 
breaking, insulator phase; at intermediate E a chiral symmetry breaking but conducting 
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phase; at large E a conducting and chirally symmetric phase. That the presence of infintes- 
simal electric field does not allow density even with a very large chemical potential and 
stops the restoration of chiral symmetry is rather surprising. The conclusion is certainly 
correct within the analysis that we have performed. It is possible though that additional 
stringy physics should be present near the IR boundary in this limit to explain the sink 
for the energy the E field is injecting. Such physics could potentially change the phase 
structure at low T. 

Finally we have explored the full E, fi,T volume at fixed B to show how these phases 
are linked. The phase diagram is summarized in Fig 1(d) with the transition boundaries 
marked. 

The variety of phase structure and transition type in such a simple theory is remarkable. 
Whether these results can serve as an exemplar for other gauge theories either qualitatively 
or quantitatively remains to be seen but they certainly suggest a rich structure of phases 
will be present in many gauge theories. 
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